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Abstract 

We address the possibility of taking advantage of high accuracy gravitational space 
experiments in the Solar System and complementary cosmological tests to distin- 
guish between the usual general relativistic theory from the alternative modified 
Newtonian dynamics paradigm. 
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1 Introduction 



The current observational status indicates that gravitational physics is in a 
considerabl e agreement with Einst ein's theory of General Relativity (GR) 



( [Willi . 120051 ; iBertolami et all l2006bl ); however, there are some theoretical and 
experimental issues that shed doubt on its status of the ultimate description 
of gravity. 



Most theoretical difficulties are related to the strong gravitational field regime, 
together with the ensuing spacetime singularities. It is believed that a possi- 
ble way of overcoming these troublesome questions lies in the quantization of 
gravity. Unfortunately, the success of modern gauge field theories in describ- 
ing the electromagnetic, weak, and strong interactions does not yet extend to 
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gravity and, in fact, the fundamental pillars of physics, Quantum Mechanics 
and General Relativity, are incompatible with each other. Furthermore, funda- 
mental theories that attempt to include gravity lead to new long-range forces. 
Even at a purely classical level, Einstein's theory does not provide the most 
general way to establish the spacetime metric, and alternative metric theories 
that have been put forward allow for violations of the Equivalence Principle, 
modification of large-scale gravitational phenomena and variation of the fun- 
damental "constants" , which are fairly constrained phenomenologically. Thus, 
the derived predictions serve well to motivate new searches for very small de- 
viations from General Relativity; these quests should clearly include further 
gravitational experiments in space, including laser astrometric measur ements 
flLATOR Collaboration fTurvshev et"ail.l2005lJTurvshev et alll2004al lb1ld\ high- 
resolution lunar laser ranging (LLR) (IMurphy et all 120021 ) and lo ng range 



track ing of spacecraft using the for mation flight concept, a s proposed (IPioneer Colaboration (Dittus et ; 
20051 ) to test the Pioneer anomaly (j Anderson et all 120021 ) ; a broad er review of 
fundamental physics experim ents in space can be found elsewhere (jBertolamil . 



20061 : iBertolami et all l2006af ). 



From an observational standpoint, recent cosmo logical observations indicate 
that, if GR is indeed correct, then one must also assume that most of the 
energy content of the Universe lies in some unknown forms of dark matter 
and dark energy that may permeate much, if not all spacetime. Indeed, re- 
cent Cosmic Microwave Background Radiation (CMBR) WMAP three year 
data( Spergel et aD . 2006 ) tells us that the Universe is well described by a flat 
Robert son- Walker metric, with an energy density of the Universe fairly close to 
its critical value, p c = 3Hq/8ttG ~ 10~ 29 g/cm 3 , where Hq ~ 71 km s^Mpc' 1 
is the Hubble expansion parameter at present. Furthermore, CMBR, Super- 
nova and large scale structure data are consistent with each other if, in the 
cosmic budget of energy, dark energy corresponds to about 73% of the critical 
density, while dark matter to about 23% and baryonic matter to only about 



As is widely known, dark matter was firstly suggested in 1933 by Zwicky, 
from a study of the observed motion of the peripheral galaxies of the Coma 
cluster; this indicated a discrepancy betwe en the total m ass, as inferred from 
its brightness or total number of galaxies (iTrimbld . 119871 ). which appeared to 
indicate that the total amount of mass in the cluster is about 400 times more 
mass than expected. This led to the postulation of a yet although unaccounted, 
non-luminous matter, which became known as "dark matter" . This hypothesis 
was also supported by the differential rotation of our galaxy, firstly discussed 
by Oort in 1927, the flatness of galactic rotation curves, large sc ale data and , 



more recently, by the observation of the so-called "bullet" cluster (j Clowe et al. 
2006h . 



Returning to current times, one finds that the amount of new data and related 
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issues coincides with recent progress in high-precision measurement technolo- 
gies for physics experiments in space, ranging from spacecraft navigation to 
time transfer, clock synchronization, weight and length standards. Hence, the 
concerned physicist is no longer restrained to mere speculation, but has a 
chance of probing crucial matters such as the nature of dark energy and dark 
matter, effects of intermediate range forces and, ultimately, the fundamen- 
tal nature of gravity. Naturally, the evolving techniques should allow for an 
ever increasing precision in measurements, with associated refinement of vi- 
able models, and exclusion of unrealistic ones. Eventually, it is expected that 
an ultimate "theory of everything" , engulfing both Quantum Mechanics and 
General Relativity will finally shed some light over the cosmological issues in- 
volving the origin and destiny of the Universe. In what follows, we discuss the 
current status of one of the most widely discussed alternatives to GR, namely 
the Modified Newtonian Dynamics and its underlying fundamental theory, the 
so-called Tensorial- Vector-Scalar theory. For the discussion of the implications 
of scalar tensor theories an d modifications designed to mimic dark en ergy the 
reader is directed to Refs. (jBertolami et all l2006bl ; iBertolamil . 120061 ) 



2 Modified Newtonian Dynamics and Tensorial- Vector-Scalar the- 
ory 



2. 1 Modified Newtonian Dynamics 



Although the current paradigm of GR endowed with dark matter and dark 
energy components is highly successful in describing observations on larger 
scales, and seeding very enticing theoretical ideas, it is important to study 
competing alternatives, in particular those which attempt to account for ob- 
servations without the dark matter component. For insta nce, a possibility 



invol v es the putative running off the gravi tational coupling (jBertolami et al 



19931 ; iBertolami and Garcia-Bellidol . Il996l ). Another, well discussed alterna- 
tive, goes by the name of Modified Newtonian Dynamics (MOND), and is 
based in a mod ification of Poisson's equation for the gravitational potential 
( IMilgroml . Il983l ). Hence, instead of the traditional form, MOND postulates 
the following equation: 



AnGp 



where p is the density of barionic matter and a$ ~ 10~ 10 m s~ 2 is the scale at 
which Poisson's equation fails; specifically, one regains the traditional form for 
accelerations a = — a , through the Milgrom function /x, which values 



3 



//(x) ~ x for and fi(x) m 1 for x ^> 1. 



This phenomenological prescription provides with an alternative solution for 
the flattening of the rotational curves of galaxies. Indeed, far away from the 
central region, the gravitational potential should be small enough so that 
MONDian behaviour takes over, and a simple algebraic computation shows 
that, instead of the inverse square law for gravity, one obtains the desired r _1 
behaviour, yielding v c (r) constant. Furthermore, a direct consequence of this 
scenario is that the luminos ity scales with y%, which is in agreement with the 
empirical Tully-Fisher law (IMilgroml . Il983l ; IZhao and Famaeyl . 120061 ) . 



Given that it was designed to solve these and other issues specifically, MOND 
is a purely phenomenological model, thus lacking a suitable support in the 
form of a fundamental theory, covariantly formulated in terms of the variation 
of a specific action functional. To avoid this caveat, several attempts to gen- 
eralize MOND were undertaken, which have so far settled in what is known 
as Tensorial- Vector-Scalar theories. 



2. 2 Tensorial- Vector-Scalar theory 



The proposed Tensorial- Vector-Scalar (TeVeS) ( Bekensteinl . 2004 ) theory as- 
sumes, asides from the metric g a p and matter fields </?j (with the associated 
energy-momentum tenor T a p, the presence of a vector field U a and two scalar 
fields, cr and (p. The vector field is everywhere timelike, with the normalization 
condition U a U a — — 1; furthermore, the vector field couples with normal mat- 
ter through a physical metric g a p = e~ 2( ^g a p — 2U a Up sinh(2(f>). The action of 
the theory is written as S = So + Sy + S$ + Sm, with 



S G 
S v - 

S s - 
Sm : 



J <r/u-s) 1/2 d 4 x 



K 
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£(^)(-s) 1/2 d 4 
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:-g) 1/2 d 4 



X 



(-g) 1/2 d 4 x , 



and K, k and / are constants specific of the theory, A is a Lagrange multiplier 
that sets the timelike normalization condition for the vector field, F is a free 
function and h af3 = g af3 — JJ a JJ° . Notice that there is no kinetic term for the 
scalar field cr, whose equation of motion provides an algebraic relation between 
its value, the free function F and the dynamic scalar field <f>. 
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There are a number of criticisms concerning the action of the theory: firstly, 
the presence of a non-dynamic o is somewhat ad hoc, although it is argued 
that the scalar action resembles that of a self-interacting complex scalar field 



in a particular limit (IBekensteinl . 12004 ) . This suggests that a is not a fun- 



damental field, but only an artificial tweak of the theory in order to enforce 
some specific behaviour. The choice of physical metric 

ga/3 ig 

also somewhat 

odd, since it is assumed that the vector field U a couples directly with mat- 
ter. Furthermore, there is no coupling between the scalar fields and matter. 
Also, asides from the dynamic and physical metrics, a third tensor h a @ also 
appears, which is somewhat similar to the physical metric, but appears only in 
the scalar action; this reveals that the vector field couples to the non-dynamic 
scalar field only, another arbitrary feature. Finally, the free function F has no 
rationale apart from the fact that it is specified in order to obtain the desired 
features. In the overall, the action for TeVeS appears as an odd concatenation 
of different tensorial, vector and scalar quantities, with no apparent inter- 
nal consistency, asides from the fact that it provides MOND with some more 
underlying supporting theory. 

The full demonstration of the connection between TeVeS and MOND is outside 



the scope of this text, and can be found elsewhere (IBekensteinl . I2004J ) . It can be 
shown that MOND behaviour arises if one assumes a spherically symmetric, 
quasistatic case, that is, one endowed with a weak potential and where only 
slow motion occurs. For this identification to be complete, Milgrom's function 
jj, is obtained from a specific choice for the free function F . This yields a rather 
complex form for the latter quantity, having an asymptote when its argument 
goes to unity, and diverging as it approaches infinity. The function F(x) also 
lacks definition when its argument varies between 1 < x < 2, thus dividing it 
into two different branches; the branch x < 1 is relevant in quasi-stationary 
systems, while x > 2 concerns cosmological issues. Given that F appears as 
a potential function for the non-dynamic field a, this increases the difficulty 
associated with the latter, as previously discussed. 



3 Testing MOND and TeVeS 



3. 1 Cosmology 



In what follows, we briefly describe the current status of observations, and 
its implications in testing if MOND and TeVeS can be a viable alternative 
to the standard paradigm of GR and dark matter. Before dwelling into its 
natural testing grounds, we quickly address some cosmological issues, namely 
the current accelerated expansion of the Universe and the initial structure 
formation. These are implemented in the context of TeVeS by assuming that 
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the vector field has only a temporal component (due to isotropy) and resorting 
to the x > 2 branch of the free function F(x). Also, ho mogeneity requires 



that the scalar field depends solely of the cosmic time t (IBekensteinl . 12004 



Skordis et all 120061 ). 



Regarding the accelerated cosmic expansion, it is found that a cosmological 
constant term can be included in TeVeS by adding an appropriate integration 
constant to the free function F; this has no effect in an astrophysical context, 
but yields the required repulsive force that count eracts the expected decel- 
eration of the rate of expansion (IBekensteinl . 120041 ) . However, its inclusion is 
purely phenomenological, and offers no explanation as to the value of this 
quantity, nor it provides any clue concerning the much debated cosmological 
constant problem. 



The issue of initial structure formation is more evolved; indeed, standard cos- 
mology assumes that structure formation occurred due to the condensation o f 
initial fluctuations, enhanced by the presence of dark matter (jTrimbld . 119871 ) . 
In the case of TeVeS, which does not assume this component to be present, 
the situation is still under study, but proponents of the theory usually rely 
on TeVeS also mani festing MOND behaviour in the early Universe scenario 
( ISkordis et all 120061 ). More specifically, instead of density fluctuations of the 
dark matter component, it is argued that fluctuations of the scalar field <fi also 
give rise to baryonic structure formation. However, thes e resu lts have been 
recently challenged by a numerical study (jPointecouteaul . 120061 ) . Concerning 
other evolutionary stages of the Universe, namely inflation, the radiation and 
matter eras, it is found that the scalar field does not interfere significantly 
with standard cosmology, as long as its magnitude does not change abruptly 
between the end of one era and the onset of the next, and is globally between 
0.0007 < < 1 (IBekensteinl . l2004h . 



3.2 Galactic Clusters 



Recent observations of galactic clusters have been interpreted as a strong indi- 
cation of the validity of the GR with dark matter paradigm. This assessment 
is based on the gravitational profile reconstructed from gravitational lensing 
data, which clearly shows two str o ng; peaks, coinciding with the centers of the 



colliding galaxies (IGavazzi et all 120031 ; IClowe et all 120061 ) . Indeed, most of 



the "visible" matter in a galaxy is in the form of plasma, with stars and dust 
accounting for only about 20% of the total; colliding plasma distributions tend 
to merge, yielding a gravitational profile with a single peak, while "normal" 
matter retains its integrity, so that its gravitational profile should show the 
individual peaks of the colliding objects. Given that the reconstituted profile 
for the two galaxy clusters is in accordance with the latter, one concludes that 
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a dark matter component dominates each galaxy's mass distribution. Since 
known matter amounts only for a small proportion, dark matter is what ac- 
counts for the gravitational profile. There are, however, some tenacious efforts 
towards an alternative interpretation of this result, trying to implement the 
MOND hypothesis in this context. This resorts to an already known feature 
of MOND (and, fundamentally, from TeVeS): that this theory does not com- 
pletely rule out the need for dark matter, but still enables a reduction of its 
proportion by a factor of five (with the remainde r possibly bein g made of 



neutrinos or some unaccounted "normal" matter) (ISandersl . Il998l ). Hence, it 
is suggested that MOND can also account for the twin peaks gravitational 
profile , by assuming the p resence of neutrinos with a mass of approximately 



2 eV (lAngus et all . l2006al lbh. Two criticisms can be put forward: firstly, the 



assumed neutrino mass is very close to the known allowed upper bound of 
0.07 eV < m u < 2.2 eV, so its validity must be reassessed as tighter results 
become available. Secondly, and possibly more crucially, the reconstruction of 
the observed gravitational profile based on MOND resorts to various linear 
combinations of "MONDian" potentials <pi . This is somewhat troublesome, 
given that the inherent non-linearities of the underlying TeVeS theory may 
indicate that, even if each of the galaxies allows for a MONDian behaviour 
away from its central region, one cannot treat their collision as a simple super- 
position of individual potentials. Moreover, the chosen potentials are not fully 
accounted for, with no clear justification given for the values of the fitting 
parameters or the algebraic form of the potentials themselves. Also, uncer- 
tainties on the mass distribution, gravitational lensing statistics and spatial 
symmetry of galaxy mergers do not allow for a clear relation between the fit- 
ting para meter giving the st rength of the scalar field 0, and the typical value 



k ~ 0.03 (iBekensteinl . [200J), which appears in the free function F of TeVeS 
and sets the order of magnitude of this scalar field in a spherical symmetric 
case, through 



kGM 
4irr ' 



where <p c is the cosmological value of 0, to which <fi(r) goes asymptotically, as 
the metric goes to its Minkowski form. Furthermore, a thorough analysis of 
galaxy merging in the context of MOND should eventually resort to TeVeS 
as a more fundamental approach to tackle the complexity of the spatial su- 
perposition of two matter distributions. A mere fit of parameters based on 
plausible, but somewhat arbitrary potentials, offers only a superficial solution 
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3.3 Solar System 



It is our opinion that implications of MOND in the Solar System should be 
better examined. Most probably, this lack of interest is due to the preconcep- 
tion that the condition V</> <C ao does not occur within the vicinity of the 
Sun (as can be checked by direct computation), and, as a result, MOND-like 
behaviour should be nonexistent or highly attenuated. However, there are two 
avenues of research that may be pursued: the first concerns the study of regions 
where the above condition may be realistically found, that is, by assuming not 
only a simple central body problem, but also take the remaining planets and 
objects into account. Specifically, one can deal with a two-body problem, and 
analyze the region close to the equilibrium point V4> = 0; a partial numerical 
study of this issue shows that there is a MONDian "bubble" of non-negligible 
size at such location, wh ich could provide interesting clu es as to the validity 



of the MOND paradigm (IBekenstein and Magueijd . 120061 ) . Although direct ef- 



fects are below current measurement capabilities, it is found that an object 
located within this bubble would experience an anomalous acceleration which, 
in principle, could be measured in the near future. 

In parallel, one can assume a simpler central body problem and, while as- 
suming that no MONDian regime occurs in its vicinity, perturbative effects 
should arise in the standard Scharzschild metric. This is clearly su itable for an 



analysis within the framework of the PPN formalism (jWilll . l200ll ). In short, a 
simplified version of this formalism states that, in a isotropic reference frame, 
the expansion of the metric to second order in the gravitational potential reads 

2GM /2GM\ 2 2GM 

(7oo = -1 + P( , <fe = l + 7 • 4 

r V r / r 



The PPN parameters (3 and 7 are related to fundamental physical properties 
of the theory under scrutiny: the latter quantifies the spatial curvature pro- 
duced by unit mass, while the former indicates the amount of non-linearity in 
the superposition law for gravity. These can be used to compare and distin- 
guish GR from competing theories of gravitation, with the former displaying 
the reference values /3 — 7 = 1. 

It was first thought that such a procedure would yield results which made 
MOND indistinguishable from GR, that is, that the PPN parameters (3 and 



7 do not shift from unity (IBekensteinl . I2004T ) . However, this assumed that the 



vector field U a had only a non- vanishing temporal component which, by virtue 
of its normalization condition, tracks the metric element g 00 and, after some 
algebraic work, leads only to a redefinition of Newton's gravitational constant 
G; this prescription for the vector field is allowed by its equation of motion , 



and argued as a plausible choice for its form. However, a later work (IGiannios . 
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20051 ) showed that one can also obtain another solution for the vector equation 
of motion, where U a has also a non-vanishing radial coordinate, related to the 
temporal one through the normalization condition. Contrary to the previous 
null deviations from the GR values for the PPN parameters, this yields a case 
which differs from standard GR0: 

= 1 + ^ + X + < M 3 + ^7±\/^7 + 5 l , 7 = 1 ■ (5) 



It should be pointed out that, in our opinion, the result (3 ^ 1, asides from 
its clear theoretical interest, is more physically motivated that the previ- 
ous (3 = 7 = 1 case. This is so because the latter is based on the Ansatz 
U a = (U°, 0,0,0), while the former takes U a = (U°, U\ 0, 0): while both are 
consistent with the vector field equation of motion, the purely temporal case 
seems somewhat artificial. Indeed, although the velocity-field of known mat- 
ter admits only a temporal component, the scalar field <p, which pervades all 
space, has a velocity-field with both temporal and radial components. Hence, 
a sensible assumption seems to be that the vector field also has temporal and 
radial componen ts. Moreover, if indee d one has (3^1, the superposition prin- 
ciple assumed in Angus et al. ( 2006al lbl) for the summation of potentials is not 



exact, which further casts doubts concerning the ability of MOND to account 
for the gravitational profiles of galactic clusters. 

Recall that the current experimental data shows an impressive agreement with 
General Relativity, with the most stringent bounds 

(3 - 1 = (1.2 ± 1.1) x 10" 4 , 7-1 = (2.1 ±2.5) x 10" 5 , (6) 



arisin g respectively from the Cassini's 2003 radiometric experiment (IBertotti et al. , 



20031 ) and from limits on the Strong Equivalence violation parameter, 77 



4/3 — 7 — 3, tha t are found to be r\ = (4.4 ± 4.5) x 10 4 , as inferred from LLR 



measurements (IWilliams et al.l . I2004J ) . 



Hence, the assumed value of k ~ 0.03 is immediately excluded, which poses an 
issue concerning the choice of the free function F, since this quantity specifies 
the order of magnitude of that parameter of TeVeS. Also, one obtains a weak 
bound on K, which must be smaller than 10 -3 ; unfortunately, there are no 
other tight constraints on this quantity to be measured against with, since the 
latter is related to the strength of the vector action in the full TeVeS action, 
and the vector field U a does not play a dynamical role in cosmology or galactic 
structure. In any case, it is clear that more stringent bounds of j3 may render 
relevant constraints for TeVeS. 



1 We have recently and independently obtained a very similar result, unaware that 
it had been derived earlier. 
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4 Discussion and Conclusions 



In this study we have discussed the current status of the Modified Newtonian 
Dynamics alternative to General Relativity, based on the Tensorial- Vector- 
Scalar theory; we address both its theoretical foundations, as well as recent 
observations that question its applicability, attempting to cast a circumstanced 
criticism of both. We conclude that, although MOND still remains a logical 
possibility, it is at the verge of being excluded from the observations, unless 
one is willing to exploit some of its theoretical oddities, most of which unac- 
counted for. Hence, although it may still be too early for its utter dismissal as 
a viable alternative to the paradigm of GR, we believe that the case for MOND 
has grown increasingly weaker, favouring the more canonical approach that 
assumes the presence of a dark matter component in the Universe. Hopefully, 
the future will hold further experimental results and theoretical advances in 
both competing theories, leading to a more profound understanding of gravity. 
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